Effects of Anisotropy in QED3 from Dyson- Sch winger equations in a box 
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We investigate the effect of anisotropies in the fermion velocities of 2+1 dimensional QED on the critical 
number A^^ of fermions for dynamical mass generation. Our framework are the Dyson-Schwinger equations 
for the gauge boson and fermion propagators formulated in a finite volume. In contrast to previous Dyson- 
Schwinger studies we do not rely on an expansion in small anisotropies but keep the full velocity dependence of 
fermion equations intact. As result we find sizable variations of A^^ away from the isotropic point in agreement 
with other approaches. We discuss the relevance of our findings for models of high- Tc superconductors. 
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I. INTRODUCTION 

Quantum electrodynamics in two space and one time di- 
mensions (QED3) has been studied intensely in the past. 
QED3 is a comparably simple quantum field theory which 
nevertheless displays a rich structure of nonperturbative phe- 
nomena such as (a variant of) confinement and dynamical 
chiral symmetry breaking plus the associated generation of a 
fermion mass. In condensed matter physics, QED3 and vari- 
ants thereof have been suggested as effective theories to de- 
scribe strongly interacting systems such as high temperature 
superconductorsia. or grapheno^. 

In the context of high-temperature superconductivity the 
chirally broken and symmetric phases of QED 3 are associated 
with the antiferromagnetic spin-density wave state (SDW) and 
the pseudogap phase of the cuprate compound described as 
'algebraic Fermi liquid'ii^. Here, the fermions of QED3 rep- 
resent chargeless quasi-particles, which are low-energy ex- 
citations close to the nodes of the (d-wave symmetric) gap- 
function. These spinons interact with collective, topological 
excitations of the gap that can be described by a U{1) gauge 
theory. Furthermore, since the motion of the quasi-particles is 
confined to the two-dimensional copper-oxide planes in these 
systems one ends up with quantum electrodynamics in (2+1) 
dimensions. 

There is a long history of debate on the critical number Nj 
of fermions, where the isotropic version of QED 3 undergoes 
the transition from the chirally broken into the chirally sym- 
metric phase. At zero temperature, if N'j > 2 then underdop- 
ing drives the superconductor directly into the SDW phase, 
while if N'j < 2 the system will first go through the pseu- 
dogap phase. Indeed, in the continuum studies mostly using 
Dyson-Schwinger equations''^ this number has been found 
in the range of A'^^ = 3.5 — 4 using Ward-identity improved 
fermion-photon vertices'*. This range seems to include the 
critical number extracted from lattice gauge theory—^, pro- 
vided the (large) corrections due to the finite lattice volume 
have been taken into accoun t '^''^ . 

When comparing with the experimental situation, however, 
isotropic QED3 may not be the appropriate choice. In fact 
there are sizable anisotropies in the cuprate compounds due 
to large differences in the Fermi velocity vp and a second 
velocity vd related to the amplitude of the superconduct- 
ing order parameter. This fermionic anisotropy can be as 



large asii X = vp /v/\ ^ 10. Furthermore, both veloci- 
ties are also different from the speed of light, c^, which is 
the third scale in the theory. Thus in the effective theory 
all three Euclidean axes have to be distinguished and it is 
an important question how and whether the ciitical number 
of fermion flavors Nji depends on {cs,vp ,va)- In the lit- 
erature, this situation has begun to be explored in a num- 
ber of approaches. First studies using renormalisation group 
and Dyson-Schwinger equations'-- used an expansion only 
valid for small anisotropies. There it has been concluded that 
anisotropy is an irrelevant phenomenon in the renormalisation 
group sense and does not lead to drastic changes in N^. In 
contrast, lattice gauge theor y '^''^ and a recent study of the 
strength of the fermion-photon couplingi^ provided evidence 
for the opposite picture; according to these works anisotropy 
is relevant and its increase leads to a decrease of the critical 
number of fermion flavors A^^. This opens up the possibility 
that Nji may go down from Nji = 3.5 — 4 in the isotropic case 
to A^^ < 2 for values of anisotropy relevant for the high-Tc 
superconducting systems. However, further studies are neces- 
sary to pin down quantitative predictions for this effect. 

In this work we reexamine the situation of anisotropic 
QED3 from the perspective of Dyson-Schwinger equations. 
We go beyond the small anisotropy expansions used in 
Refsjiii^ and solve the Dyson-Schwinger equations for the 
nonperturbative, anisotropic fermion propagator numerically 
using on the one hand the l/Af/-approximation and on the 
other hand an improved approximation of the photon propaga- 
tor as input. These technical details are summarized in section 
mi Our results, presented in section [Till confirm the findings 
from the lattic o'^''^ and from the fermion-photon coupling 
strengtbi^. We discuss qualitative aspects of the dependence 
of Nf on the anisotropies and identify meaningful approxima- 



tion schemes for the DSEs. 
section HVl 



We conclude with a summary in 



II. TECHNICAL DETAILS 

A. The Dyson-Schwinger equations in anisotropic QED3 

The details of the formulation of anisotropic QED3 have 
been discussed in Refs. fllfTsllTil Here we only summarize the 
most important elements using the notation of Ref. HI The 
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fermions are represented in the four component spinor repre - 
sentation obeying the CUfford algebra {7^, 7^} = 2 (J^j,. We 
consider Nf fermion flavors and the anisotropic fermionic ve- 
locities vp and VA which will be included in the metric-like 
factor gi.pi/. In the high-temperature superconducting (HTS) 
system, the d-wave symmetry of the gap function on the Fermi 
surface leads to four nodes, which can be grouped into two 
distinct pairs^. This is accounted for by the index i = 1,2. 
The "metric" is given by: 
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VF^ 

va'- 



(1) 



and 
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va' 



(2) 



vp- 



It enters the Lagrangian in the form of: 



J, I'M 



(<9p + i a^) 



(3) 



As in isotropic QED3, we can define chiral symmetry using 
73 and 75. With massless fermions, the Lagrangian there- 
fore has a U{2Nf) "chiral" symmetry, which is broken to 
SU{Nf) X SU{Nf) X [/(I) X U{1) if the fermions become 
massive. The order parameter for this symmetry breaking is 
the chiral condensate which can be determined via the trace of 
the fermion propagator 5_F,j (p) = / ^|^e*^^*(i'j (f)'!'^ (0)). 
The corresponding expression for the photon propagator is 
given by D^,{p) = J ^e'P-{A^{x)A,{0)). The Dyson- 
Schwinger equations for the fermion and photon propagators 
are diagrammatically represented as shown in Fig. |l(a)| and 
Fig. |l(b)| In Euclidean space- time, they are explicitly given 
by 



SI 



i=l,2 



(27r) 



rTr 



(5) 



with the momentum k defined by the difference p — q and i = 
1,2. Here, the renormalization constant Zi of the fermion- 
boson vertex, ( p, g ), is included. 

As we consider an anisotropic space-time for fermions, we 
have to take care of dressing each Dirac component separately 
and also include the node index. In order to keep track more 




FIG. 1: The diagrammatic representation of the Dyson-Schwinger 
equations for the fermion (a) and gauge boson propagator (b) . Wig- 
gly lines denote photon propagators, straight hnes fermion propaga- 
tors. A blob denotes a dressed propagator or vertex, whereas a dot 
stands for a bare fermion-photon vertex. 



easily of the general structure of the equations, we introduce 
some shorthands denoted by 

Pi — -Pm 9i Pu (6) 

and 

p^^i :— A^^i{p)p^, (no summation convention !), (7) 

Where A^^i denotes the vectorial fermionic dressing function 
at node i. 

The anisotropic expressions for the dressed fermion and 
photon propagators then read 

Splip) = S,(p)+i VsI'^T^Pm.- (8) 

D^. (p)" = p2 (^6^. - ^) + n^. (p ) , (9) 

where Bi denotes the scalar fermion dressing function at node 
/ and n^,y ( p ) the vacuum polarization of the gauge boson 
field. As the free gauge bosons do not involve fermionic ve- 
locities, the bare propagator still is the isotropic one. 

In order to close the Dyson-Schwinger equations for the 
propagators we need to specify the fermion photon vertex F*^. 
This vertex satisfies its own Dyson-Schwinger equation in- 
volving a fermion four-point function of which very little is 
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known to date. Strategies to explore the fermion-photon in- 
teraction in the isotropic case therefore centered around ver- 
tex constructions based on higher order perturbative expan- 
sions, multiplicative renormalizability and the Abelian Ward 
identity One of the main results of Ref. |8i, however, is 
the stability of the (isotropic) critical number of fermion fla- 
vors with respect to the dressing of the vertex: the difference 
between the values for a bare vertex and much more sophisti- 
cated constructions is smaller than five percent. This finding 
somewhat justifies the use of simple vertex constructions also 
in the case of anisotropic space-time, where the complexity 
of the DSEs is large and any simplifications are highly wel- 
come. In this work we therefore employ a vertex construction 
which only takes into account the leading term of the Ball- 
Chiu vertex^- (IBC-vertex), generalized to anisotropic space- 
time. It reads 



(10) 



where no summation convention is used and p, q are the 
fermion and anti-fermion momenta at the vertex. This vertex 
is known to be a reasonable approximation to the full Ball- 
Chiu vertex^, which solves the Ward-identity exactly. We will 
use the vertex Eq. ( fTOb in the quark DSE and also compare 
results with simpler truncation schemes such as the 1 /Nf ap- 
proximation or the Pisarski schemes^. 

In order to make the photon DSE tractable we will use 
two different approximations: the approximation and a 

more sophisticated model based on the results of Ref. HI In the 
limit of isotropic QED3 there has been extensive work in the 
l/Nf approximation^'^, which has been found to agree qual- 
itatively with more elaborate approaches'^. In the anisotropic 
case the leading order l/Nf approximation has been pursued 
in Refs. fUfTsi where in addition an expansion in terms of small 
anisotropics has been performed. This is not the case here: we 
keep the full dependence of the equations on all three veloci- 
ties {cs,vp ,va)- As a result, one may expect to get a rough, 
but qualitatively correct picture of the effects of anisotropy 
close to the transition from the chirally symmetric to the chi- 
rally broken phase. 

The vacuum polarization in leading order large- A^/ expan- 
sion is given diagrammatically in Fig.|2l The integration can 
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FIG. 2: The diagrammatic representation of the Dyson-Schwinger 
equation for the vacuumpolarization of gauge bosons in leading order 
large-AT/ expansion. 



be performed analytically and one obtainsiii^: 



n.(p) 



A^f 



(11) 



Despite its merit of simplicity, there is a serious drawback 
of the leading order 1/A^/ -expansion, which makes it desir- 
able to go beyond: it leaves the values of the fermion vector 
dressing functions Ai = 1 at all momenta. In the isotropic 
case, this behavior is in marked contrast to the one of more so- 
phisticated approximations. There, infrared power laws in A, 
the photon polarization 11 and the quark-photon vertex appear, 
which are characteristic for the chirally symmetric phased. In 
addition, in the anisotropic case, non-trivial values of the A- 
functions at zero momenta serve to define renormalized ve- 
locities (cf^, VF ^, VA It is not clear whether this is pos- 
sible within the 1/A^y -expansion. In order to go beyond this 
approximation we therefore utilize the results of Ref. § in 
the isotropic case and in the symmetric phase the photon vac- 
uum polarization is well approximated by the following ex- 
pression: 



1 



I^1+2k J^2 g2 



(12) 



This expression includes the 1 /A^/ -limit for k = 0. In general, 
however, this anomalous dimension is non-zero and depen- 
dent on the details of the vertex truncation. In Ref. [U the value 
of K for the vertex truncation Eq. (flOb at the critical number 
of fermion flavors A^^ has been determined as k = 0.0358. In 
principle, this number may depend on the fermion velocities 
in the anisotropic case. Nevertheless, due to our current in- 
abiUty to calculate this dependence we will keep k constant in 
all calculations. We therefore arrive at the anisotropic gener- 
alization of (fT2l i 



f 



1 



1 



16 vpVA 



m 



(13) 

with pf defined in Eq. (|6]l. 

In general, we wish to emphasize, that great care has to 
be taken to use consistent approximations in the photon and 
fermion DSEs. Either one has to use the 1/A^/ approxi- 
mation in both DSEs, leading to the well known isotropic 
result^ A^^ 3.24 , or one has to use the vertex Eq. (fTOb in 
both DSEs and solve everything self-consistentlyS', leading to 



A^' 



/ 



3.56. Below we will show that we reproduce this last 
result approximately using our model Eq. (fT3] l for the pho- 
ton. It is, however, inconsistent and therefore wrong to use 
the l/A^/-approximation for the photon polarization Eq. (fTTI) 
and the vertex Eq. ( fTOb in the quark-DSE. This choice leads 
00 and therefore prohibits any chiral restoration. 



to A^' 



/ 

We also wish to emphasize that to our mind, meaningful re- 
sults cannot be extracted from models that violate translation 
invariance of the photon propagatori2i2^. 
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Having specified our approximation (flST l of the photon, we 
now proceed with the fermion-DSE. In order to extract the 
scalar and vector dressing function of the fermion we take ap- 



propriate traces of the fermion-DSE and rearrange the equa- 
tions to end up with: 



B^{p) 



= Zo 



Z.e' 



P/. J (27r)" 



B^ {qy + iq-^r 

I 



(14) 



(15) 



where again we have no summation for the external index /i. 
Furthermore we have used the Ward-identity Zi = Z-i be- 
tween the fermion wave function renormalization factor Z-2 
and the one for the vertex, Z\. The form of the ''metric^ leads 
to several symmetries that will simplify the evaluation of the 
DSEs. In general, we find the relations: 



Bl{po,Pl,P2) = B2{po,P2,Pl) , 
{P0,Pl,P2) = A^,^2{P0,P2,Pl) , 



n" {P0,P1,P2) 
11^°{P0,PUP2) 



(pO,P2,Pl) , 
n20(po,P2,Pl). 



(16) 
(17) 
(18) 
(19) 



One then ends end up with eight anisotropic equations that 
need to be solved self-consistently, in comparison to only 
three equations for an isotropic space-time. This increase 
in the number of equations results from the component-wise 
dressing of the fermion momentum and from the anisotropic 
vacuum polarization of the gauge bosons. In our case using 
the l/A^/-approximation or our model ([13] ) in the photon sec- 
tor we remain with four equations. These reduce to two in the 
isotropic limit. 



B. The DSEs on a torus 

Although our approximations of the photon serve to sim- 
plify the system of DSEs dramatically, the numerical proce- 
dure to solve the remaining equations for _Bi, Aq.i, ^11,^2.1 
is still considerable due to the presence of the anisotropy. In 
particular, the usual transformation to hyper-spherical coor- 
dinates and exploitation of symmetries therein is of course 
no longer feasible. One possibility to circumvent this prob- 
lem is to change the base manifold of the system from three- 
dimensional Euclidean space to that of a three-torus. In 
Lorenz-covariant systems such a treatment has been explored 
extensively in the context of QCD^^ '^^ and also for QED3A2.. 
Here we generalize this treatment to the anisotropic case. In 
fact, the formulation of DSEs on a torus is ideally suited 
for this endeavor, since a formulation in Cartesian instead of 
hyper-spherical coordinates is most natural. 

On a compact manifold, the photon and fermion fields have 
to obey appropriate boundary conditions in the time direc- 
tion. These have to be periodic for the photon fields and 



anti-periodic for the fermions. For computational reasons it 
is highly advantageous, though not necessary, to choose the 
same conditions in the spatial directions. We choose the box 
to be of equal length in all directions, Li = L2 = L3 = L, 
and denote the corresponding volume V ^ L^. Together with 
the boundary conditions this leads to a discretization of mo- 
mentum space. Thus all momentum integrals appearing in the 
Dyson-Schwinger equations are replaced by sums over Mat- 
subara modes. For the fermions this amounts to 



d^q 
(27r)3 



(...) 



1 

l3 



E (•••). 



(20) 



Tii,n2,n3 



counting momenta q^^ = si'^'^/ + l/2)ei, where 

ii are Cartesian unit vectors in Euclidean momentum space. 
For the photon with periodic boundary conditions the momen- 
tum counting goes like qn = 3(27r/i)(ni)ei. The re- 
sulting system of equations can be solved numerically along 
the lines described in Ref . [l3l 

Of course, working on a torus with a finite volume intro- 
duces artifacts. These have been studied in some detail in 
Ref.dH and found to be sizable, see Fig. 7 of this reference. 
In addition, nonlinear volume effects have been found for very 
large volumes in Ref. Q The reason for these large effects is 
the large separation of scales in QED3. Besides the natural 
scale = 1 of the coupling with dimension one, there is 
the scale of dynamical mass generation which even well in 
the broken phase is at least two orders of magnitude below 
the natural scale. This is a situation totally different to, e.g., 
QCD where both scales are of the same order In the iso- 
topic theory, the large volume effects in QED3 explain why 
lattice simulations^^ find critical numbers of fermion flavors 
of the order of Nr « 1.5, while continuum calculations in the 



DSE framework^ find N 



f 



3.5. The DSE calculations on 



the torus explained these two results as a consequence of the 
finite lattice volumd^^. We will briefly address the correspond- 
ing volume effects in the anisotropic case in section llllCI We 
find indications, that the size of the volume effects is hardly 
affected by changes of the fermion velocities, at least in the 
range relevant for our study. Thus while the absolute number 
for Nj on the torus needs to be corrected, its variation with 
(cs, vf , WA ) is nevertheless meaningful. We therefore adopt 
a similar strategy as the authors of Ref. [Ill and concentrate on 
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the variations rather than the absolute value of NJ. 



III. NUMERICAL RESULTS 

In this section we will discuss our numerical results for 
three different approximation schemes. For the fermion DSE 
these are (i) the Pisarski scheme""', (ii) the leading order l/Nf 
expansion and (iii) our most elaborate scheme using the vertex 
Eq. ( fTOb and the photon model Eq. ( fTSl ). Before we present our 
results for the isotropic and anisotropic cases at finite volume, 
we briefly discuss the corresponding infinite volume results in 
the isotropic case. In general, besides the chiral condensate 
another suitable order parameter for chiral symmetry break- 
ing is the value of the scalar fermion dressing function at any 
point in momentum space; here we use B{p) at the smallest 
momentum available. 



A. The critical Nj^ in the isotropic case 

In the Pisarski scheme, analysed in detail in Ref. |23l the 
fermion-photon vertex is set to the bare one and the scalar 
dressing function of the fermion is set to a constant, B{p) — 
B{Q) = m. This rough approximation scheme is well suited 
to analytically demonstrate that chiral symmetry is broken in 
QED3. As a result one obtains 



5(0) 



cae 



(21) 



where c is a positive constant and a = Nf. Clearly, chi- 
ral symmetry is broken and remains broken for any value of 
Nf. Thus N'j — > 00, which is in contrast to the findings in 
more elaborate approximation schemes. Therefore, although 
the Pisarski scheme served as a nice tool for demonstrating 
broken chiral symmetry in the first place, it is not well suited 
to extract information on chiral symmetry restoration. We will 
see below that this remains true in the anisotropic case. 

Let us now focus on the other two approximation schemes 
(ii) and (iii). In the infinite volume limit the dependence of 
5(0) on Nf has been determined in Ref. together with re- 
sults for more elaborate vertex truncations. Here we com- 
bine the result for our photon model Eq. (fTSl ) together with 
the vertex approximation Eq. (fTol ). In Fig. [3] we display our 
results. For all three truncations we clearly find an exponen- 
tial decrease of the order parameter with Nf, a situation de- 
scribed as Miransky-scaling-*^ in the literature. (Note that uni- 
versal power law corrections to Miransky scaling have been 
discussed in Ref. and Refs. therein. For QED3 the quan- 
titative aspects of these corrections have not yet been stud- 
ied.) The critical number Nj has been determined analytically 
for the and the self-consistent IBC-vertex truncations: 



V 



3.24 and N 



c,lBC 



value for our model TV^ ™"' 



del 



3.56. The corresponding 
is a little lower than the one for 
the 1/A^j^-approximation, as can be inferred from Fig. [3] 

The reason why 7v^'™°'*<^' ^ N^'^^^ is the following: the 
photon in the self-consistent IBC-approximation is a con- 
stant in the infrared for Nf < Nj and develops an infrared 
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FIG. 3: The order parameter for chiral symmetry breaking, B{0), 
plotted as a function of the number of fermion flavours Nf in the 
isotropic case 1 = Cs — vf — va ■ Shown are results for the self- 
consistent leading order -approximation and the IBC-vertex, 
Eq. diet , compared with the results for our model Eq. ( 113) for the 
photon together with the vertex Eq. dlO) in the fermion-DSE. 



power law for Nf > Np In our model, the power law is al- 
ready present for Nf < N'j, which results in missing infrared 



strength and therefore explains why TV 



c^model 
f 



< N 



c,lBC 
f 



This is, however, only a small quantitative effect which is ir- 
relevant for the qualitative study presented in this work. 



B. The critical in the anisotropic case 

In this section we discuss the general anisotropic case with 
Cs = I and {vp , va ) varied. We first investigate the Pisarski 
approximation, which has been used in Ref. to argue for a 
couple of general points. Again, similar to the isotropic case, 
we can determine the behavior of the order parameter ana- 
lytically, at least for vp = va ■ The details are reported in 
appendix|A] We find 

m = B(0) = caexp (-TT^Nf^/vp^ - l) / (22) 
(2y/vF^ - 1 + 2(2 + vp^) arctaii(VwF^~T) 

which is different to the corresponding expression reported in 
Refill Our result has the correct isotropic limit and it gen- 
eralises the finding from the isotropic case: in the Pisarski 
approximation (and for d = 3 dimensions) chiral symmetry 
is always broken, i.e. Nj — > 00. We therefore find this ap- 
proximation not to be suited in general to extract quantitative 
information on the variation of Nj with {cs,vp ,va)- 

We now proceed with our two more elaborate approxima- 
tions. In Figs. |4] and 12] we show our results on a torus with 

= 1, L = 600/e^ and 20'^ momentum points. The values 
of N^ are represented by contour lines and a coiTesponding 
color code. Note that the 'edges' in the contour lines are due 
to the sparseness of our grid in some regions; in general we ex- 
pect smooth lines in the vp -va -space. For vp — va = 1 we 
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FIG. 4: Critical number of fermion flavours A'^^ determined in the 
-approximation and plotted as a function of the fermion veloc- 
ities {vf ,va)- 




FIG. 5: Critical number of fermion flavours Nj determined from 
the photon model Eq. J13b and the fermion-photon vertex Eq. JlOt 
plotted as a function of the fermion velocities {vf ,va)- 



can read off the corresponding values for Nj in the isotropic 

■del 



case. We find N'j'^^^'^ 



1.0 and Nf"""" 



1.1. Clearly, 

these values are affected by large finite volume effects, as al- 
ready discussed above. When we compare these with the cor- 
responding values in the continuum calculation we find that 
the latter are roughly a factor of three larger This has to be 
kept in mind when discussing Figs. |4] and |5] 

For non-trivial values of vp and wa we find sizable devia- 
tions of Nj from the isotropic case. We also observe that the 
two truncation schemes behave very differently. In general the 
l/A^/-approximation shows an increase of the critical num- 
ber of fermions with respect to deviations from the isotropic 
case in all directions. This finding is in agreement with early 
expectations'^ from Dyson-Schwinger equations, performed 



in the l/iV/-expansion and furthermore in an expansion of 
small anisotropics. Thus what we have shown here is that the 
results of Ref. [Ts] survive also for large anisotropics. There 
remains, however, the question of the reliability of the 1/Nf- 
expansion. We have already discussed above, that the leading 
order l/A^/-expansion contains the approximation Ai ~ 1, 
which is certainly not valid in general. Our modified approx- 
imation scheme Eq. ( fTOb and Eq. (fTSl ) avoids this problem. 
Indeed, our results for in this model show a different be- 
havior. Here we find a maximum of for small values of 
(f i? , tja ) at or around {vp ,v^) = (0.4, 0.4) with decreasing 
behavior for smaller and larger velocities. In the region where 
both velocities are smaller than c,, = 1 the critical number of 
fermions never goes to zero. This is different for large veloci- 
ties: in fact if one of the velocities is larger than a certain crit- 
ical value of the order of vp ,va ~ 2, the system is always 
in the chirally symmetric phase. Certainly the precise num- 
ber, where this happens, is expected to be severely affected by 
volume effects. Nevertheless, this qualitative behavior is in 
good a gree ment with the corresponding evaluation performed 
in Ref. llsl, where the criterion of the strength of one-photon 
exchange between fermions was used. Our results therefore 
corroborate their findings. (It is perhaps amusing to note that 
in Ref. llSlthe statement was made that it is 'almost impossi- 
ble' to analyze mass generation in the presence of anisotropics 
using Dyson-Schwinger equations. Naturally, in the light of 
this work we do not agree with this statement.) As argued in 
detail in Ref. [T^ the picture indicated by Fig. |5]is also sup- 
ported by lattice calculation o'^i'^ . In the light of this agree- 
ment we tend to the statement that the 1/A^j -approximation 
picture of Fig.|4]is misleading and has to be discarded. 



C. Volume effects and the case vf /va 
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Let us now briefly come back to the volume effects asso- 
ciated with our calculation on the torus. In Ref. (TU these ef- 
fects have been analyzed in the isotropic case and tori as large 
as = 512'^ have been used. In the anisotropic case the 
numerical effort to determine is much larger than in the 
isotropic case and the N = 20'^ calculations reported above 
are already demanding in terms of CPU time. A detailed and 
concise study of the volume effects using large tori is there- 
fore beyond the scope of the present work. Nevertheless we 
wish to demonstrate, that the assumption used above, namely 
that the size of the volume effects do not too strongly depend 
on the size of the anisotropy, may be correct. To this end we 
show the results for four selected points on the vp ^ va- 
axis for N = 30'^ in Fig.|6l Clearly, the larger volume results 
in a larger value for Nji, as visible from the plot. Compar- 
ing the volume effects at the isotropic point with the value at 
Vp = 1.5 we find quantitative differences, but certainly not 
too different to invalidate our assumption from the last section 
on a qualitative level. Finally we also observe, that the line of 
demarcation from 7^ to = moves strongly with 
volume. Clearly, further studies are necessary to pin down the 
size of the volume effects to such an extent, that extrapola- 
tions to the infinite volume limit along the lines of Ref. [ij are 
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calculation in Fig.|2l The filled circles in the plots are our nu- 
merical results, whereas the black Une is a linear fit to the data 
points. Clearly, the dependence of Nj on ua is well approx- 
imated by this fit. Whether this linear dependence remains 
intact in the infinite volume hmit should be studied in future 
work. 



0.5- 



IV. SUMMARY AND CONCLUSIONS 



0.0- 



1.0 



1.5 



2.0 



FIG. 6: Four selected points on the vp = va -axis determined with 
two different tori with A'' — 20^ and A'^ = 30^ corresponding to box 



lengths of Le^ = 600 and Le^ 



900. 




FIG. 7: Selected results for the critical number of fermion flavors at 
{vf ,va) = (lOuA , '^A ) from a torus with Le^ = 600. The error 
bars correspond to our step size when searching for Nf. 



In this work we have proposed an approximation for the 
Dyson-Schwinger equations of QED3, which is well suited to 
investigate the effects of anisotropics on the critical number of 
fermion flavors for chiral symmetry restoration. This problem 
is important for the description of high Tc cuprate supercon- 
ductors using QED3 as an effective theory. Our results support 
corresponding ones from lattice calculation a'^ '^ and contin- 
uum studies of the fermion-photon interaction strength^. In 
our best approximation scheme we find a maximum of the 
critical number of fermion flavors around vp = t'A ~ 
0.5. From this point, Nj is decreasing in all directions in the 
{vp ,Va )-plane. Further studies, however, are necessary to 
pin down the quantitative aspects of this behavior Our re- 
sults are still preliminary in at least three aspects. They are 
affected by large volume effects, which have been studied in 
detail only at the isotropic velocity point'^. Our calculation is 
certainly affected by truncation errors in the fermion-photon 
vertex, although we have argued from comparison with more 
systematic calculations at the isotropic point**, that these ef- 
fects are on a moderate level of the order of twenty to thirty 
percent. Finally, our treatment of the anisotropic photon DSE 
remains preliminary and a more self-consistent calculation is 
certainly desirable. Nevertheless, we believe that the quali- 
tative aspects of our results are correct and provide a further 
step forward to understand the connection between QED3 and 
experimental studies of high Tc superconductors. 



possible. To forego the problem of finite volume effects, one 
may also attempt to solve the equations directly in the infinite 
volume continuum limit. We will report on first results in this 
direction elsewhere. 

When comparing with experiment, values of {vp , ua ) of 
the order of X = vp /va ~ 10 may be of particular inter- 
est; Ref . [Til reports A = 14 for optimally doped YBa2Cu307 
and A = 19 for Bi2Sr2CaCu208 compounds. In addition, it 
is important to know the relation of {vp ,va) to the speed 
Cs related with the gauge boson. Here, we recall that Cg is 
not the speed of light in the cuprate compound, but the speed 
characteristic for the vortex-antivortex excitations that are rep- 
resented by the gauge field in the effective theory. It is thus 
possible and even anticipated that vp > Cg, as discussed in 
Ref. [m To our knowledge, the precise relation of vp and 
Cs is not yet determined from experiment. We therefore show 
results for (IOwa , va ) in a typical range accessible in our 
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Appendix A: Anisotropic Nf in the Pisarski approximation 

In this appendix we give some details to the dependence 
of the fermion mass on the number of flavors in the Pisarski 
approximation. For vp = ua , we insert the large-iV/ vac- 
uumpolarization, Eq. (fTTl i. into the scalar fermion dressing 
function given implicitly in Eq. For equal fermionic ve- 
locities, we make use of the inherent symmetry and switch to 
cylindrical coordinates. With the abbreviations 
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A = 2vF ^z'^i^/r^ + z^a^ + 16z^{A^/r^ + + a)) + 8r^{8^/r^ +z^ + a + VF "^{S^/r^ + + 3a)) 

B = r^VF^z^{{l + VF^)^\/r'^ + z^a^ + &z^{S\/r^ + + '24:VF^\/r^ + z^ + Idvp'^Vr^ + z^ + {2 + 5{vf^ + VF'^))a)) 

C = 2r^{vF V?-^ + z^a^ + 4:Z^{8\/r^ + z^ + a + VF^il6\/r^ + z^ + 24wi. + z^ + {5 + 5vf ^ + vf ■*)«))) 

V = Nf{r^ + z^){r^ +VF'^z'^)^{8\/r^T^ + a) 

£ = 64^"^ + r2(64(l + VF ^)z^ + 8^/ + z'^a + vf '^a{8^/r^ + + a)) + vp '^z'^{64:Z^ + a{16\/r^ + z^ + a)) 



we amve at 



27r oo oo 



B{p) 



(27r)3 



difdrdz 



16 B{p) raiA + B + C) 



(Al) 

At this point we perform Pisarski's approximation^ and set 
the scalar fermion dressing function B{p) equal to a constant 



J 



r 



B{p) = m for all momenta in the range a ^ p ^ m. Other 
momentum contributions are neglected in the integration. Fur- 
thermore, since r < a and q < awe only keep terms with the 
largest power of a in the numerator and denominator of the 
integrand. This leads to the more compact equation 



1 7. ?. } . 16mr{2r'^ + r^{l+VFYz^ + 2vF^z^) 



(27r)3 



dip dr dz 



Nf{r'^ + z2)2 (r2 + VF ^z^)^ 



m m 

As a next step, we perform the z integral and again neglect the terms not proportional to a. We end up with 

27r a 



1 



m 



dip / dr 



8m (Vi'F ^ — la + {2 + a) arctan[-\/wF ^ — 1]) 



(27r)3 J J Njra^JvF'^- 1 

rn 



(A2) 



(A3) 



We further perform the r and ip integral and use one last time that a is the largest scale in the problem. We finally arrive at 

2m [y/l + VF^ + (2 + VF ^) arctan(VT~rJ7/2)) 



7r2iV/Vl + VF ^ 



(A4) 



which leads to the result given in the main text, Eq. ( |22] |: 



m = S(0) = ca exp y^-Ti'^ N f^/ vf ^ - ij / [2^/vf ^ - 1 + 2(2 + vf ^) arctan( V^'f ^ - 1) 
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